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Abstract 



We study the off-diagonal Goldberger-Treiman relation (ODGTR) and 
its discrepancy (ODGTD) in the N, A, tt sector through 0(p 2 ) using heavy 
baryon chiral perturbation theory. To this order, the ODGTD and axial vec- 
tor N to A transition radius are determined solely by low energy constants. 
Loop corrections appear at 0(p 4 ). For low-energy constants of natural size, 
the ODGTD would represent a ~ 2% correction to the ODGTR. We discuss 
the implications of the ODGTR and ODGTD for lattice and quark model 
calculations of the transition form factors and for parity-violating electroex- 
citation of the A. 

PACS Indices: 11.30.Rd, 13.75.Gx 
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I. INTRODUCTION 



The Goldberger-Treiman relation (GTR) |TJJ plays an important role in theoretical 
hadronic and nuclear physics. It relates hadronic matrix elements of the weak axial current 
(the nucleon axial charge, gA, and the pion decay constant, F n ) to quantities governed by 
the strong interaction (the pion nucleon strong coupling constant, g n NN, and nucleon mass, 
m N ): 

_ 9A m N m 
9wNN — ^ UJ 

The GTR represents an approximation, since g n NN is determined experimentally at the 
point q 2 = m 2 while gA is measured close to the point q 2 = 0. In the chiral limit, the GTR 
would be exact, while in the physical world, it holds to an astonishing level of accuracy. 
The small difference between the physical value of g n NN and RHS of Eq. (P when physical 
values of gA, F m and tun are used is called the Goldberger-Treiman discrepancy (GTD). 
Physically, the GTD is driven by the explicit chiral symmetry breaking introduced by the 
non-zero current quark mass. 

Many theoretical discussions of this chiral symmetry-breaking effect have appeared in the 
literature 0-f|. Recently the GTD in the context of SU(3)lxSU(3)# chiral symmetry has 
been analyzed by Goity et al. |7|] within the framework of heavy baryon chiral perturbation 
theory (HB;\PT) ||||. These authors found that chiral loop corrections appear at 0(p A ). 
The dominant contribution comes from the low energy counterterm appearing in the 0(p 3 ) 
Lagrangian. Their result is consistent with more conventional approaches where current 
quark mass plays an explicit role 

In this work we analyze the off-diagonal Goldberger-Treiman relation (ODGTR) and its 
discrepancy for the SU(2) ir, N, A sector. As we show below, both the magnitude of, and 
theoretical uncertainty in, the off-diagonal Goldberger-Treiman discrepancy (ODGTD) is 
~ m^/A^ ~ 0.01, where A x = 47rF n ~ 1 GeV is the scale of chiral symmetry breaking. 
Consequently, the ODGTR provides a useful benchmark for both experimental and theo- 
retical studies of the axial vector N —>■ A transition form factors. In principle, the ODGTR 
can be tested using charged current reactions, such as neutrino excitation of the A, or weak 
neutral current processes, such as parity-violating (PV) electroexcitation. These processes 
are sensitive to axial vector transition form factors, which can be related to the strong 7rAA 
coupling via the ODGTR. The values for these form factors obtained from charged current 
scattering are fairly uncertain. A measurement of the PV asymmetry for neutral current 
electroexcitation will be performed at the Jefferson Lab by the GO Collaboration |TI| in 
hopes of providing a more precise determination of the axial transition form factors. The 
ODGTR also provides a check on lattice QCD and hadron model calculations of the axial 
transition form factors. From either perspective, the theoretical analysis of the ODGTR 
using HBxPT appears to be a timely endeavor. 



II. NOTATIONS 

We follow the standard HB^PT formalism |||| and introduce the following notations: 
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S = e , £ = , 7T = ~7T a T a (2) 

with F v = 92.4 MeV being the pion decay constant. The chiral vector and axial vector 
currents are given by 

= ^ + <V 



tr, r 



/r = ^re f ± e^re 

X = 2B (s + ip) 

X± = ^ t ±ex t e , (3) 

where s,p, a M , are the scalar, pseudoscalar, pseudovector, and vector external sources with 
p = p l r % and a M = a^y. 

For the A, we use the isospurion formalism, treating the A field T l { vector spinor 



in both spin and isospin space 11] with the constraint r l T^(x) = 0. The components of this 
field are 




(4) 

The field also satisfies the constraints for the ordinary Schwinger-Rarita spin-| field, 

y-r; = o and p^t; = o . (5) 

We eventually convert to the heavy baryon expansion, in which case the latter constraint 
becomes v^T't = with being the heavy baryon velocity. 

In order to obtain proper chiral counting for the nucleon, we employ the conventional 
heavy baryon expansion, and in order to consistently include the A we follow the small scale 



expansion developed in JlTJ. In this approach external energy and momenta, the A and 
nucleon mass difference 5 = hia — and 1/mjv are all treated as 0(e) in chiral power 
counting. 

The leading order HB^PT Lagrangian reads: 
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= N[iv ■ D + 2g A S ■ A]N - iTj*[iv ■ - 5 ij 5 + g x S ■ A^T^ 
+ g» NA [T?^N + Nupt\ + ?Tr[D^ED M Et + X ^ + X^} + • • • (6) 

where is the Pauli-Lubanski spin operator and uu l = Tr (r 1 ^). 

At subleading order we collect only the 7riVA interaction pieces which are relevant in the 
following discussion. 

4 1} = ^mkv^ - ^-J^N + h.c. + • ■ ■ (7) 

where 

^ = Tr^[V„A u ]) . (8) 



III. OFF-DIAGONAL GOLDBERGER-TREIMAN RELATION AND ITS 

DISCREPANCY 

It is convenient to introduce the ttNA form factor GVjva via the effective Lagrangian: 

ArJVA = -^^A WiV + h.C. (9) 

In terms of the couplings appearing in Eq. (|J) , one has 

<J-7rAfA — p > l-LUJ 

where g^ivA is the renormalized 7rNA coupling constant. We also express the matrix elements 



of the axial current between A + and proton in terms of the Adler form factors [|T~2^p~4 
< A+(p')\Al\P{p) >= A+VHQV)^ + C -^q,q, 

+ [^ + ^' A ](gAV-^)}ti(p) , (ii) 



where we have displayed only matrix elements of the neutral component, A^ = qj^s-^q, 

~J -A 4~ i i i i t -/ i ill,-, /]: iiinn.ui i- , \\ \ 



for brevity. Experimentally, one expects contributions from C§ to give the dominant effect. 



For future reference, we also define the off-diagonal charge radius, r\: 

d 



ri = 6^1nC 5 V)| g2=0 (12) 
To arrive at the ODGTR, it is useful first to contract Eq. (|TT1) with g M , yielding 



< A + (p'WAl\p(p) >= *A+V)[^V) + ^P-q 2 ]qMp) • (is) 



m 



N 



We compute the same matrix element from the amplitudes of Fig. 1. The pion pole contribu- 
tion (Fig. lb) depends on G n NA(q 2 ) and P(q 2 ), the coupling of the pseudoscalar current to 
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pions. At lowest order, one has P{q 2 ) = m 2 F n . We parameterize the non-pole contributions 
(Fig. la) in terms of a function C(q 2 ). We thus obtain 



< A+tf)\d»Al\P(p) >= -^/fiA+V) x 2 D{q l\ . q u u(p) (14) 

q z — m z + le 



with 



D{q2) = G^tfPvtf) + 2 _ ^ 2 



Equating (113) and (0), using Eq. fll5]) , and taking the limit q — > 0, leads to 



ty(o) = - < /;[- g ^g' (0) +c(o)i . tie) 



We emphasize that Eq. fll6|) involves no approximation. However, neither G^ata^) nor 
P 7r (0) is experimentally accessible. To the extent that these quantities vary gently between 
q 2 = m 2 and q 2 = we may replace them in Eq. (16) with their values at q 2 = m 2 . 



Assuming pion pole dominance and neglecting C(0) would then lead to the ODGTR. The 
off-diagonal Goldberger-Treiman discrepancy (ODGTD), A^, embodies the corrections to 
these approximations. Including A„- we have the corrected ODGTR: 



3 m N rril 

where, to leading order in light-quark masses, we have 

d 
dq 



C*(0) = _ _2 (1 - ( 17 ) 



K = ml—\nD(q 2 )\ q 2 =m 2 . (18) 



An analogous expression for the diagonal GTD case was first derived in Ref. 0. Indeed, 
our treatment here largely follows the outline of that work. 

In order to obtain A„., one requires the g 2 -dependence of both G n N&(q 2 ) and P^q 2 ) as 
well as the non-pole amplitude C(0). To that end, we first observe that since P(q 2 ) = m^F^ 
at lowest order, Cg (0) starts off as O(p ). The non-pole term C(0) generates an 0(p 2 ) 
correction, as we discuss in the following section. In principle, since P(q 2 ) is 0(p 2 ) at 
leading order, one might expect its g 2 -dependence to arise at 0(p A ). However, there exist no 
operators in the 0(p 4 ) Lagrangian ( see Ref. [0 ) which contribute to this g 2 -dependence, 
nor do the corresponding loop graphs contribute at this order. 

The g 2 -dependence of G n NA(q 2 ) requires more care. As we show explicitly below, loop 
contributions to this g 2 -dependence arise first at 0(p 4 ), and thus, for our analysis, may be 
neglected. However, in the nonrelativistic theory obtained via the heavy baryon expansion, 
the 6 3 + b 8 terms contribute to the g 2 -dependence via the factor 

y.^^-^-f . (19) 

H 2m N V ; 

Note that this term is nominally 0(p) in the small scale expansion, since m\—rn 2 N /2mN ~ 5. 
However, it contains an 0(p 2 ) contribution (the q 2 term) as a consequence of kinematics. 
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Since we derive expressions below valid in the nonrelativistic theory, we should include this 
contribution to GWa^ 2 )- 

To complete analysis of G n N/\(q 2 ), we observe that loop corrections renormalize the bare 
ttNA coupling g® N & — > g n NA at 0{p 2 ). However, the g 2 -dependence of the vertex due to loop 
corrections appear 0(p 4 ). Since we truncate at 0(p 2 ), these corrections can be neglected, 
and all we need to do is to replace g® N & by g n NA- A similar situation holds for the diagonal 
GTD, as shown in the analysis of Ref. [0. In our case this observation directly leads to the 
conclusion that the A^ and r\ are solely determined by the counterterms. 

It is useful to examine the g 2 -dependence of loop effects in some detail. To that end, we 
first classify the various diagrams contributing to the ODGTR. Diagrams (a), (e), (g), (i), 
(j) and (k) contribute to the tensor structure g^ while the remaining diagrams contribute 
to the structure q^q u - The first diagram (a) in FIG. 1 is the tree level one. The second 
diagram (b) is the pion pole contribution. Diagram (c) and (d) renormalize P^q 2 ) and 
their contribution is of 0(p 4 ) as explained above. The loops in diagrams (e) and (f) contain 
no g 2 -dependence. Diagrams (g)-(n) are similar to each other, so we take diagram (g) as 
example. The amplitude reads 



iM, 



(a) 



g\9nNA f d D k k^S ■ qS ■ k 



F 2 J (2tt) d k 2 - m 2 n + ze 

X ^k+ k2 ^v-(k + q) + ^^^ (20) 

where q is the external momentum and we include the leading recoil correction in the nucleon 
propagator. According to HBxPT, the recoil corrections may be included perturbatively, so 
we expand the baryon propagators in (pH) as follows: 



iM {9) 



d u k k^S-qS-k 1 1 
(2ir) D k 2 — m 2 + iev-kv-k + 5 



r k 2 — (v ■ k) 2 (v-q) 2 — 2k-q v ■ q, 

The first term inside the square brackets generates a g 2 -independent contribution of 
0(p 2 ). Upon integration, the terms in the integrand containing explicit factors of q generate 
an additional factor of v ■ q/m^ relative to the leading term. According to Eq. (|19D , this 
factor contains a g 2 -dependent term which goes as —q 2 /2m 2 N . Thus, the g 2 -dependence of 
this integral occurs at 0(p 4 ). Similar arguments hold for the other loops in diagrams (h)-(n). 



IV. THE LOW ENERGY COUNTER TERMS 

Consider first A,,-. We collect the 0(p 3 ) low energy counterterms which may contribute 
to A,,: 

41 = -T^Tt[D„ X _]*N + ^T;[D U , r\N (22) 

x x 

+^|^ 75 [ X _^] J iV + h.c. + --- 
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where [D^x-Y = Tr{ yf-D^, X-}} e ^c The ellipses denotes other 0(p 3 ) terms which do 
not contribute to A n . Detailed expressions of these terms can be found in Ref. |16|]. After 
carrying out the heavy baryon expansion, the third term in Eq. (^) is of 0(p 5 ), where one 
power of p arises from a factor of p/rriN generated by the 275 tensor structrure. Also the 
third term contains two pion fields. So its contribution to A^ involves with one additional 
loop and is further suppressed by 1/A X . In other words, this piece can be neglected. 

Since we obtained our general expression for A^ using matrix elements of d^A^, we may 

(3) 

deduce its dependence on the q by varying C C j. with respect to the pseudoscalar source, p\ 
To that end, we use the chiral Ward identity of QCD 



T 



d^lqi^^q] = m[qi^ 5 r l q\ 



with rh = m " + md . Moreover, 



(23) 



5C 



gvy S T q 



QCD 



5p % 



From Eqs. 



, g4| ) and the leading-order relation x % - = 4iB op 1 we obtain 

5£hb x pt 



Equations (^3| , P^ , P^D and (|2"5D then imply that 



C£{q 2 ) + 



C£(q 2 



m% 



-q 



''N 

+{h + 6 8 ) 




r 

2 \9* 



JVA 



m, 



m 



N 



2mArA, 



x 



(24) 



(25) 



(26) 



where we have used IB^m = m%. With Eq. ( 18) we arrive at the off-diagonal GTD to 
0(p 2 ): 



A, 




2ci ~h + b 8 f A, 



girNA 2(yf 7r ATA \TnN 



(27) 



The ODGTD - whose scale is of order (m n /A x ) 2 ~ 0.01 - depends on three low-energy 
constants: g^NA, c\i and 63 + b 8 (we count the latter as a single constant). Since we have 
scaled out explicit factors of 1/A X in C C ' T , we expect these constants to be order of order 
unity. In fact determinations of g n NA and 63 + b$ from ttN scattering in the resonance region 
yield m 

g nNA = 0.98 ± 0.05 
63 + b 8 = 0.59 ±0.10 . 

Were C\ also to be of order unity, we would expect A^ to be of order a few percent. This 
magnitude for A n is consistent with previous estimates |5|]17| . As in the diagonal GTR the 
ODGTR should hold to within a few percent accuracy, as a consequence of chiral symmetry. 
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Consider now the leading g 2 -dependence of C§(q 2 ). Since loops do not contribute to 
the g 2 -dependence of C§(q 2 ) at 0(p 2 ) we need consider only the tree- level contributions 

(3) (3) 

generated by C c ^. They are most easily obtained by considering the dependence of L CT on 
the pseudovector source a 1 - 



A' 



5C 



HB X PT 



(28) 



We then arrive at 



v 9ttNA + (&3 + b 8 



m, 



2 1 1 



2m N A x 



1 2 
q 

x x 



(29) 



so that 



^x 



C-2 



+ 



63 + b 8 / A 



9nNA \m N 



(30) 



where we have dropped higher order contributions (e.g., corrections of order 5/mjv)- From 
Eq. (p6|) we also conclude that 



gm^ATAL^z^ - 6r i] + °(Q 



2 ml' 



^m^F^G^NAi-^: 



6r\] + 0(q 



2 mV 



(31) 



Note the low-g 2 behavior of the induced off-diagonal pseudoscalar form factor is completely 
determined (once r\ is known) , since it is expressed in terms of the physical and measurable 
parameters as can be seen from the second line in Eq. fl3~T|). 



V. IMPLICATIONS FOR EXPERIMENT AND THEORY 

In principle, an experimental test of the ODGTR could be carried out by drawing upon 
precise measurements of C5 (0) and G^na^D- A value for C5 (0) has been obtained from 
charged current neutrino scattering from hydrogen and deuterium ||18|| : 



C^(0) = -±=(2.0±0.4) , (32) 

where the prefactor is due to relative normalization of charged and neutral current ampli- 
tudes. 

For the strong 7rNA form factor, one may rely on the analysis of ttN scattering given in 
Ref. |I6"|1 , which gives 



G^ NA (ml) = 11.6 ±1.3 . (33) 

Substituting this result into Eq. (|17]) and dropping the correction A n yields the leading-order 
ODGTR prediction for C£(0): 
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C 5 A (0) Lo . = 0.93 ±0.10 



(34) 



A comparison of this value with the experimental result in Eq. (|32|) leads to an experimental 
constraint on the ODGTD: 

A cxp = _ Q 2A ± Q 3 ^ ( 35 ) 

where the error is dominated by the experimental error in (7^(0) . 

Alternately, one may draw upon the older analysis of the K-matrix for pion photopro- 
duction |19|J20[1 in the A resonance region to obtain 



G« NA (ml) = 14.3 ±0.2 , (36) 

which implies 

A^ p = 0.01 ± 0.2 . (37) 

In both cases, the value of A| xp is consistent with zero and, thus, in line with our 
expectations that the ODGTD be of order a few percent at most. At present, however, the 
uncertainty A^. xp is an order of magnitude larger than one would like in order to test this 
theoretical expectation. Since this uncertainty is dominated by the error in C£(0), it would 
be advantageous to reduce this uncertainty through more precise form factor measurements. 

Such measurements could also reduce the present uncertainty in r\, which has been 
determined from charged current neutrino scattering data. An empirical parametrization of 



C*5 (q ) obtained from this data gives [^1 



1 + 1.21-^-2-3 — 

C 5 V) = C*(0) - 2G a f- q2 (38) 



V 1 M\> 

with Ma = 1.14 — > 1.28 GeV. From this parameterization, one would deduce 

,1.21 2 



f = (— + w ) = (1-82 -> 2.14)GeV- 2 (39) 



Accordingly we determine 



c 2 = -(3.1 -> 3.5) (40) 



While the value for c 2 is consistent with expectations that it be of order unity, its uncertainty 
is roughly 10%. 



Parity- violating (PV) electroexcitation of the A, as approved to run at Jefferson Lab |T0 
will provide new, precise measurements of the axial vector — > A amplitude at a variety 
of q 2 points. At first glance, this program of measurements could yield a determination of 
both C^(0) and r\. However, the extraction of these quantities from experiment requires 
resolution of two theoretical issues. The first involves the overall normalization of the axial 
vector amplitude and, thus, the determination of C^(0). The normalization - which could be 
obtained from a fit to the measured g 2 -dependence |22| - is strongly affected by electroweak 
radiative corrections, R^, as discussed in detail in Ref. [23]|. As emphasized in that work, 
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these corrections are theoretically uncertain, as a result of nonperturbative QCD effects, 
and the corresponding uncertainty could be on the order of 10-20% relative to the tree-level 
amplitude. The radiative corrections always come in tandem with axial vector amplitude 
for PV electroexcitation and cannot be determined independently (e.g., by proper choice 
of kinematics or target). Thus, they introduce an intrinsic, theoretical uncertainty in the 
extraction of C5 (0) from this process. Given the estimated size of the uncertainty, it appears 
unlikely that PV electroexcitation will improve upon the result in Eq. fl52]). 

Nevertheless, determining the normalization of the axial vector amplitude via the Jeffer- 
son Lab measurement would be interesting from another perspective. Because the theoretical 
uncertainty in the ODGTD is considerably smaller than both the current experimental error 
in C^(0) as well as the estimated theoretical uncertainty in R^, one might use the ODGTR 
prediction for C^(0), in tandem with the normalization of the axial vector amplitude ex- 
tracted from PV electroexcitation, to determine R^. Recently, the study of axial vector 
electroweak corrections has taken on added interest in light of the results of the SAMPLE 
experiment | 25| , which imply that the magnitude of Ra for elastic, PV electron scattering 



may be considerably larger than implied by theory [26]. Understanding these corrections 



could have important implications for the interpretation of other precision electroweak mea- 
surements, such as neutron /3-decay |2"7]] , so it would be of interest to study them in both 
the elastic and inelastic channels. 

A second interpretation issue involves the g 2 -dependence of the PV asymmetry and, 
thus, the determination of r\. In contrast to the situation for elastic, PV electron scattering 
- where the PV asymmetry vanishes linearly with q 2 at low-|g 2 |, the asymmetry for PV 
electroexcitation contains a g 2 -independent term. In the framework of Ref. [23], this term 
is characterized by a low-energy constant d&. On the scale of the expected asymmetry, the 
magnitude of the d& contribution could be significant, particularly at low-|g 2 | where one 
would want to determine r\. In order to determine the latter reliably, one also requires 
knowledge of d&. 

The second issue could, in principle, be resolved through a measurement of A 7 , the asym- 
metry for PV photoproduction of the A. Since v4 7 is proportional to d&, and since chiral 
corrections to the asymmetry are small, its measurement could remove the a^-related uncer- 
tainty in PV electroexcitation. Thus, measurements of both Ay and the PV electroexcitation 
asymmetry at a variety of q 2 points could yield values for r\, d&, and R^. 

New, precise neutrino scattering experiments would complement this program. Since 
neutrino scattering probes of the axial vector transition amplitude are free from the large and 
theoretically uncertain radiative corrections entering PV electroexcitation, such experiments 
could, in principle, provide a theoretically clean determination of C^(0). 

Finally, we observe that the ODGTR could provide a theoretical self-consistency check 
on lattice QCD and hadron model computations of the axial vector iV — > A transition form 
factors. While there exist lattice calculations of the electromagnetic N —>■ A amplitudes, 
the axial vector amplitudes remain to be computed. The lattice electromagnetic amplitudes 
appear to differ significantly from experimental values, and it would be useful to have a 
corresponding comparison in the axial vector channel. Historically, a variety of hadron 
model calculations of C^(0) have been performed, with predictions generally lying in the 
range 0.8 — > 2.0 (see Ref. [Q for a compilation). Those lying near the lower end of this 
range are most consistent with the ODGTR, based on the value of GWaI^) from Ref. 
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16fl . For example, the quark model calculation of Ref. predicts C^(Q) in terms of 
and the nucleon and A masses: 



^(O)q.m. = (-^xM 9A = 0.87 . (41) 

1.175V3 VmA + mjv/ 



The leading order ODGTR prediction is given in Eq. (|34|), where the uncertainty is domi- 
nated by the error in G n N&(ml) obtained from Ref. HT6| . Thus, the quark model appears to 



be consistent with the expectations derived from chiral symmetry and the latest analysis of 
strong interaction data. Having in hand similar agreement with future lattice calculations 
would be similarly satisfying. 
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FIG. 1. Relevant Feynman diagrams in the derivation of off-diagonal Goldberger-Treiman 
relation and its discrepancy. The filled circle denotes the pseudoscalar or pseudovector source. 
The double, solid and dashed lines correspond to the delta, nucleon and pion respectively. 
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